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ABSTRACT 

The prediction of collisions amongst N rigid objects may be re¬ 
duced to a series of computations of the time to first contact for all 
pairs of objects. Simple enclosing bounds and hierarchical par¬ 
titions of the space-time domain are often used to avoid testing 
object-pairs that clearly will not collide. When the remaining pairs 
involve only polyhedra under straight-line translation, the exact 
computation of the collision time and of the contacts requires only 
solving for intersections between linear geometries. When a pair 
is subject to a more general relative motion, such a direct collision 
prediction calculation may be intractable. The popular brute force 
collision detection strategy of executing the motion for a series of 
small time steps and of checking for static interferences after each 
step is often computationally prohibitive. We propose instead a less 
expensive collision prediction strategy, where we approximate the 
relative motion between pairs of objects by a sequence of screw 
motion segments, each defined by the relative position and orienta¬ 
tion of the two objects at the beginning and at the end of the seg¬ 
ment. We reduce the computation of the exact collision time and of 
the corresponding face/vertex and edge/edge collision points to the 
numeric extraction of the roots of simple univariate analytic func¬ 
tions. Furthermore, we propose a series of simple rejection tests, 
which exploit the particularity of the screw motion to immediately 
decide that some objects do not collide or to speed-up the predic¬ 
tion of collisions by about 30%, avoiding on average 3/4 of the 
root-finding queries even when the object actually collide. 
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Figure 1: Collisions between a moving an object A, the block 
with a through hole, and a static object B, the L-shaped extru¬ 
sion. A moves from right to left along a screw and its intermedi¬ 
ate instance shows it at the moment of first collision with B for 
two different configurations: a face of A collides with a vertex 
of B (top) and an edge of A collides with an edge of B (bottom) 

Introduction 

The calculation of collisions amongst moving 3D objects and be¬ 
tween moving objects and static obstacles has challenged anima¬ 
tion experts and medical engineers for more than two decades [1, 
7, 10, 38, 27, 32]. It has also been extensively studied in robotics 
[36], We focus here on rigid bodies, and more precisely on poly¬ 
hedra, and do not address collisions of deformable models [48, 2]. 
When complex 3D motions are involved, the poses (position and 
orientation) of each moving object are usually evaluated using a se¬ 
ries of small time increments. At each stage of this simulation, the 
transformed instance of each object is tested against the instances 
of other objects using a static interference test [37]. Simple contain¬ 
ing bounds [44], convex decompositions [3, 18, 19], hierarchical 
models [4, 9, 42, 22, 25, 26, 29, 35, 39, 41], and minimal distance 
computation [6, 15, 23, 37, 43] and tracking [13, 14, 5, 24] tech¬ 
niques have been used to reduce the frequency and the complexity 
of the interference tests [17, 20, 21, 31]. When an interference is 
detected, the last interval may be refined adaptively, until an accu¬ 
rate interval around the initial collision time is isolated. Even when 
instantaneous velocities or bounds on velocities are used to esti¬ 
mate the duration of collision-free intervals [11, 21], the number 
of interference test that are necessary to guarantee that all colli¬ 
sions are detected imposes limits on the performance of these col¬ 
lision detection approaches. By contrast to the above approaches, 
which are based on series of static interference tests, we propose a 
collision prediction approach, in which we compute the time and 
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location of collisions directly from the relative motion of pairs of 
objects. Instead of computing swept volumes[36, 28, 34, 46, 49], 
or intersecting four-dimensional models swept by the moving ob¬ 
jects in the space-time domain [12], we detect all occurrences of 
face/vertex and edge/edge collisions [16], and report the first one 
to occur. When both objects move along straight-line translations 
[10, 11] or both are rotating around the same axis [47], these direct 
collisions calculations may be reduced to the evaluation of a series 
of linear or quadratic expressions. 

For more general motions where linear translation and rotation 
about different axes are allowed, trigonometric functions in the col¬ 
lision equation cannot be removed if the rotation angle and the 
translation displacement are interpolated as a linear function of a 
single parameter. In [17, 33,47] these trigonometric functions were 
removed by nonlinear interpolation of the rotation angle, resulting 
in cubic or higher order polynomials. These nonlinear interpola¬ 
tions on angle can be considered to approximate linear ones, when 
the rotation angle is small. However, note that if the center of the 
rotation is far from the geometric center of the object, nonlinear in¬ 
terpolation on rotation produces a distorted trajectory, which may 
be quite different from the linear one. Thus, care must be taken in 
choosing the center of rotation. 

Instead, to obtain a simple formulation of the exact collision pa¬ 
rameters, we propose in this paper to approximate the relative mo¬ 
tion 1 between any two objects by a continuous series of screw mo¬ 
tion segments. For each screw motion segment, we compute the 
times of collision between all vertices of the first object and the 
faces of the other, between the faces of the first object and the ver¬ 
tices of the other, and between each edge of the first object and each 
edge of the second object. We report the smallest of these times. 
We assume of course that objects are initially disjoint. 

Furthermore, to reduce the overall computational costs, we have 
developed fast interference tests between simple three-dimensional 
or parametric bounds. We use them to quickly reject most object 
pairs, vertex/face pairs, and edge/edge pairs that do not collide. 
The collision test and the computation of the time of first collision 
for each one of the remaining vertex/face and edge/edge pairs re¬ 
quire solving a low-degree trigonometric equation in one variable. 
Once the minimal time of collision is found, we compute the colli¬ 
sion point. Our direct collision computation techniques and simple 
rejection tests could of course be combined with hierarchical ap¬ 
proaches that have been mentioned earlier. In the next section, we 
briefly explain how the screw motions are computed. Then, we 
present the simple rejection tests. Finally, we provide the details of 
the exact time to collision calculations and present implementation 
results. 

Approximating screw motions 

We propose to use a direct screw motion 2 to approximate the mo- 
tion of an object A relative to a possibly moving object B. The 
screw motion interpolates the relative pose (position and orienta¬ 
tion) of A at the beginning and ending of a given time interval, 
which we parameterize with t varying between 0 and 1. The reader 
further interested in screw theory and pose interpolation can refer 
to [28, 40, 50]. An interpolating, direct, screw motion is unique, 
except for the special case of 180° rotation between the starting 
and ending poses. By definition, it interpolates any two poses by 
a combination of a minimal angle rotation with a shortest vector 
translation. Note that the rotation axis is parallel to the translation 

^he ralative motion of an object A with respect to a moving object 
B is the motion of A in the body coordinate system of B 
2 A screw motion with samllest angle. 


vector. The interpolation is independent of the choice of coordi¬ 
nate system, thus the roles of A and B may be interchanged without 
affecting the results. 

Because the maximum error 3 of a screw approximation of a mo¬ 
tion segment is expensive to evaluate exactly for all points of an 
object, we advocate the following conservative estimation. We ex¬ 
press the discrepancy motion, D(t ), as the product of the approx¬ 
imating screw motion by the inverse of the original motion. We 
then apply D(t) to the eight vertices of an axis aligned 4 block con¬ 
taining A. The error is bounded by the maximum distance between 
the original position of each one of these vertices and its image by 
D(t), as t varies between 0 and 1. Hence, we have reduced the er¬ 
ror estimation problem to one of computing geometric bounds on 
the trajectory of several isolated points. When the error estimate 
exceeds a predetermined tolerance, the time interval is split in two 
and the motion split into two screw-motion segments continuously 
joined at the time where the two intervals meet. 

The parameters of the direct screw motion may be computed eas¬ 
ily for each segment. Assume that at instant t, the pose of object A 
is represented by a transformation that is the composition of a rota¬ 
tion Ra( 0, which is represented by a 3 x 3 matrix, with a transla¬ 
tion by a vector va ( t ). Assume that the pose of object B is similarly 
represented by a rotation R and a translation v#(f). We can ex¬ 
press the relative motion of A with respect to B as the combination 
of a rotation R (t) with a translation by a vector \(t). R (t) may 
be represented by a 3 x 3 matrix that results from the multiplica¬ 
tion of the transpose of the matrix of R s(t) by the matrix of Ra( 0- 
The vector \(t) is the difference between va( 0> transformed by the 
transpose of the matrix of R#(f), and v#(f). 

The parameters for a direct screw motion that interpolates the 
relative position of the two objects at times 0 and 1 may be simply 
derived from the matrix R r (0)R(l) and from the vector v(l) — 
v(0). These parameters define the direction s, a fixed point p, a 
translation distance d along s, and the angle b of a rotation around 
the axis parallel to s and passing through p. 



Figure 2: Screw parameters s,p, d,b 


Computation of Screw Parameters s, p, d and b 

Suppose that two poses of an object at time time t = 0 to t = 1 are 
given and we want to find a screw motion that interpolates them. 
In this section, we briefly discuss computing such screw parame¬ 
ters. Full description of computing pose interpolating screw pa- 

3 The error is defined for each point of an object at a given time 
as the distances between its images in original and approximating 
poses. 

4 Aligned to the axes of global coordinate system. 
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rameters can be found in [45]. Let R x (t),R y (t),R z (t) represent 
respectively the first, second, and third columns of R (t) and let 
AR*^ = (R^ j3 , jZ ( 1) — R^y^O)). Fig. 3 illustrates these vectors. 
Then the rotation axis and rotation angle can be computed as 


s 

h 


s/|s| 
2 sin 


, where s = AR X x ARy + AR V x AR Z + AR Z x AR X 
-i |AR*|/2 
|sxR,(0)| 


( 1 ) 


We have expressed s as the sum of three cross-products, so as to 
guarantee that the formula works in all situations and gives the 
maximum accuracy. Note however that at most one of AR X , AR } ,, 
or AR, can be null and therefore at most two of the cross-products 
can be null. Therefore, if AR X and AR-y are both not null, we could 
also compute s as AR X x AR-y. 

After computing the screw axis direction vector s and rotation 
angle b , a point on the screw axis p and translation d can be com¬ 
puted. Consider a point o in the object. Let o (t) be the position of 
o at time t. Since poses of time t = 0 and t = 1 are given, o(0) and 
o(l) are also given. From Fig. 3, d and p can be simply computed 
by 


d = (o(l) — o(0)) s 
P = \ fo(l)+o(0) + 


sx (o(l)-o(O)) 
tan (b/ 2) 


( 2 ) 


Sphere bounding swept 



Figure 4: Trivial rejections: cylinder/sphere collision 



Rejection Tests for Entire Objects 

• We first check the collision between the bounding sphere, 

Sph(oB,rs ), of B and the sphere, Sph(os,rs ), that encloses 
the volume swept by a bounding sphere of A. The symbols 
o and r with the appropriate subscripts represent respectively 
the centers and radii of these spheres. As is illustrated in 
Fig. 4, center and radius of Sph(og , r$ ) are computed as = 

p' + d/2 s and rs = + \Jd 1 j 4+ |oa — p'| 2 . The rejection 

condition is \o s — o#| > r# + rs 

• If the two bounding spheres are intersecting, we check whether 
the sphere around B and the infinite cylinder centered around 
the axis of the screw and containing the Sph(oA , r a) intersect. 
First, the distance between the screw axis and the center of 
the bounding sphere of B is computed. The rejection condi- 


Figure 3: Computation of screw parameters s,p, d,b 


Early Rejection Tests 

In this section, we propose a series of simple tests for identifying 
situations where collision is clearly impossible. We assume that, 
in a direct pose interpolating screw motions, the magnitude of the 
rotation angle b is less than 7T, without loss of generality. When 
expressed in a relative coordinate of B, object A moves along a 
relative screw motion and object B is static. We assume for sim¬ 
plicity that the boundaries of both objects have been triangulated 
and that the objects are initially disjoint. (The latter conditions 
may be tested using a static interference test.) We must test for 
vertex/triangle, triangle/vertex, and edge/edge collisions. Since the 
triangle/vertex case can be converted to a vertex/triangle case by 
swapping the roles of A and B, or equivalently by inverting s, there 
are only two different cases to consider. If the two objects have 
n and m vertices, 0(nm ) pairs of entities need to be tested. We 
describe here quick rejection tests. They could be combined with 
hierarchy based algorithms, such as those for example based on 
sphere hierarchies [30]. 


tion is 

|o^p-o^p-ss| > |o A -p'|+r A + r B (3) 

If |p' — Oa | is larger than r^, the following condition can also 
be used as a rejection test, identifying the cases where B lies 
inside a shrunken version of the cylinder. 

|o^p-o^p-ss| < \o A -p'\-r A -r B (4) 

• Finally, the planes that cap the finite cylinder that contains 
the volume swept by the sphere around A are considered. If 
| (og — os) • s| >J/2 + rA+rg, there is no intersection. 

Vertex/Triangle rejection tests 

If the above tests fail to reject the pair of objects, we test for colli¬ 
sion between the individual pairs of elements of their boundaries. 
First we test all vertex/triangle and triangle/vertex pairs. As de¬ 
picted in Fig. 5, the following rejection situations for a vertex and 
triangle pair are considered. 

• After projecting the vertices of a triangle to the screw axis, 
the distances to the cylinder that contains the helix along 
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Figure 5: Trivial rejections: vertex/triangle 


which q moves are tested. If (qo — v/) • s > 0,Vz = 1,2,3 
or (qo — Vj) • s < — d, Vz = 1,2,3, then there is no collision. 

• Rejection is also possible when the triangle lies inside the 
cylinder. If | (v* — p) • s| < | (qo — p) • s|, Vz = 1,2,3, then there 
is no collision. 

• Trivial rejection is also possible when the triangle is out of 
the angle swept by the vertex as projected on a plane per¬ 
pendicular to the screw axis. From the plane equations that 
include the screw axis, qo and qi, if 0 < b < 7t, (qo — p) x 
s • (v* — p) < 0, Vz = 1,2,3 or (qi — p) x s • (v* — p) < 0, Vz = 
1,2,3, there is no collision. If n < b < In , the six conditions 
must be satisfied at the same time, i.e., if (qo — p) x s • (v; — 
p) < 0 and (qi — p) x s • (v,- — p) > 0, Vz = 1,2,3, there is no 
collision. 

Edge/edge rejection tests 

Finally, we develop rejection tests for edge/edge pairs. Assume 


• When the minimum distance between the edge bib 2 and the 
screw axis is larger than r max , there is no collision. Let n = 
s x 4/|s x lfc| be the common normal of the screw axis and 
bib 2 - Then the minimum distance is y for some t and j8 that 
can satisfy the equation p + ts + ny = bl +1^/3. The reject 
condition is |y| = |(bi -p) xs-l fc /(nxs-l*))| > r max . 

Collision Time and Contact Calculation 

As justified above, when the trivial rejection test fail, we need to 
test for the following collision cases: vertex/triangle and edge/edge. 
For each case, the collision condition is formulated as a function of 
time and is solved numerically using a Newton’s iterative method 
with the careful choices of initial guesses. 

Screw Trajectory in Analytic Form 

Each point of an object undergoing a screw motion moves along a 
helix around the screw axis, which contains a fixed point p and has 
the direction s as illustrated in Fig. 7. Our deterministic collision 
prediction is based on the Rodrigues equations [8]. Let qo be the 

q 0 = q(0) 



rsin(tb) 



Figure 6: Trivial rejections: edge/edge 


that an edge of A having vertices sl\ and a 2 is moving along a screw 
segment. Let bi and b 2 be vertices of an edge of static object B. 
For convenience, define the two radii as r m \ n = min((ai ? 2 — p) • s) 
and r max = max((ai )2 -p) • s). 

• If the projected line segment does not intersect with the cylin¬ 
drical annulus, there is no collision. The rejection condition 
is (ai ? 2 - bi ? 2 ) • s > 0 or (ai^ - b^) • s < -d. 


• When the edge bib 2 is inside the annulus, there is no colli¬ 
sion. The rejection condition is | (bi ? 2 — p) * s| < r m i n . 


Figure 7: Rotation and projection of a point qo 

initial position and q(t) be the position at time t. Let p r be an 
orthogonal projection of qo onto the screw axis and q' be qo rotated 

about s by the angle tb. We first derive p'q'= q' — p', which is the 
rotational motion of vector from p' to qo about s by the angle tb. 
Assuming |s| = 1, from Fig. 7, 

pV= q '(t) — p 7 = rcos(tb) 5-5^ + r sin(ffe) ^ XS 

v r \j) 

= cos (tb) p r qo +sin(t&) p r qo xs 

The vector from p to q'(f) can be computed as 

pq' =p'q' +(pqo s)s 

= cos(tb) pqo +sin(tb) pqo xs +(1 — cos(f^))(pqo-s)s (O 
s0(pq*O ,*b) 

where we defined a new notation 0 (a, tb) as a vector a rotated about 

the screw axis by the angle tb. From (5), q '(t) = p+ pq /= p + 
0(pqo? tb). The screw trajectory q(t) can be obtained by adding 
the translation vector tds. 

q(t) = tds + q'(t) = p + fds + 0(pqo,f&) (7) 
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Vertex/Triangle Case 

The plane that contains a triangle is defined by a unit normal and 
a signed distance d s from the origin. Any point q in this plane equa¬ 
tion satisfies d s + q • = 0. The intersection between this plane and 

the screw trajectory can be calculated by substituting equation (7) 
into this plane equation. 

0 = di + p-n i +^(s-n i ) + ^(pqo,ffc)-ii i (8) 


We can rewrite (8) in the following form 

/(f) = co + c\t + c 2 cos (tb) + C3 sin (tb) 
= co + c\t -\-Acos(tb — 6) 

Where the coefficients are given by 5 * * * * * 


c 0 = p-n s + (p(fc) -s)(s-n s ) + d s ,ci=ds-n s 

^2 =pqo-n^ + (pqo •s)(s-n 5 ) , c 3 = (pqj) xs) • n, (11) 

4 = \/ c 2 + c 3 > 0 = tan^ 1 (c 3 ,c 2 ) 

Since /(f) does not have an analytic solution, we use Newton it¬ 
erations. Consider extremal points of /(f) where /'(f) = 0. Since 
/(f) is monotonic between these points, if the signs of at two con¬ 
secutive extremal points are identical, there is no solution between 
them. If the signs are different, there is one solution. The extremal 
points computed from the condition /'(f) = c\ — Ab sin(Z?f — 6) = 0 
are 

f = (6 + (— l) m a + m7i) /b , a = sin -1 (12) 

where m is integer. Define m\ to be the smallest m such that f >f 0 - 
Then m\ can be calculated as 

=^^ 41 ) < i3 > 


m\ = min ( 2 


— 0 — cl T bto 
2iz 


where it is assumed that b > 0 without loss of generality. Once we 
find mi, we can increase it by one to calculate the next extremal 
point and then compare the signs of /(f) at these two points. If 
the signs are different there is one solution and we need to perform 
Newton iterations starting with the reflection points as the initial 
guess, which can be computed from /"(f) = —Ab 2 cos (bt — 6) = 0 
that yields bt — 0 = nn + n/2 where n is integer. We have the 
reflection points at 


0 HlC 71 /2 

b 


(14) 


We start our root-finding process with n\, the smallest n such that 
f > 0, which is computed as 


n\ 


1 e 

2 71 


(15) 


Note that in pose interpolating screw motion, be [— 7T, n\ and there 
could be up to two reflection points in the interval and we need 
to start from these two points. Also, note that Newton iteration in 
general converges in three to four iterations. 

Whenever a new helix/plane intersection point is found, a trian¬ 
gle containnment test is performed to ensure that the vertex actually 

5 These coefficients can also be derived in the local coordinate sys¬ 

tem (i,j,k,o) of the screw, where o = p',k = s, and i is parallel to 

p'qo- The screw trajectory is q(f) = |qo|cos(fZ?)i + |qo| sin(fZ?)j + 

tdk. The screw intersects plane d s + qn 5 = 0 for values of f satis¬ 

fying 

d s + (|q 0 | cos(ffc), |qo| sin(tb),td) • n s = 0 (10) 



Figure 8: Vertex/Triangle collision. 



Figure 9: Triangle/Vertex collision. 


collides with the triangle. If the intersection lies inside the triangle, 
the search is stopped. Fig. 8 and 9 show the trajectory and collision 
points in 3D on the right and the plot of the /(f) with its two roots 
in the left. Once the solution is found, only the roots that pass the 
the triangle inside test are accepted. If more than two roots passes 
the inside test, we choose the one with smallest f since only the first 
collision is interesting. 

Edge/Edge Collision Case 

Let (ai,a 2 ) be an edge of A and (bi,b 2 ) be an edge of B. Let 
\ a = a 2 — ai and 1^ = b 2 — bi. We express a point on the edge 
of A as ai + a\ a and a point on the edge of B as bi + /31 where 
0 < a, P < 1. Let q a (f, a) be the point ai + a\ a transformed by 
the screw motion at time f. Then, the collision detection problem 
amounts to finding three variables f, a and j8 for which q«(f, oc) = 
bi + Ph- Qa can be calculated by substituting ai + a\ a into (7), 

q a (t,cc) = q a (t) + a<j>(l a ,tb) (16) 

where q fl (f) and 0(1 a ,tb) are defined as 

qa(0 =v+tds+<f(m,tb) 

0(1 a ,tb) = cos(tb)\ a + sin(tb)\ a x s + (1 — cos(tb))\ a s s 

Note that 0(1 a ,tb) is the transformed \ a at time f. The collision 
condition is 


qa{t,cc)=q a {t) + a<j>(l a ,tb) = b + pi b (18) 

To eliminate a , we take the cross product with 0(1 a ,tb), yielding 
x 0(l fl , tb) = b x 0(1 a ,tb) 4- Ph x 0(l fl , tb) and then we take the 
dot product with 1^ resulting in q a x 0 (\a,tb) • 1^ = b x 0 (l fl , tb) • 1^. 
Using the vector identity uxv-w = u- vxw, this can be written as 

/(f) = (q fl -hi) • (0(lfl,ffc) x\ b )=0 (19) 
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Figure 10: Edge/Edge collision. 


Using vector identities 6 and assuming |s| = 1, f{t ) can be simpli¬ 
fied to the following form 

f{t ) = co + ( c\ + C3 1) cos (tb) + (C2 + C41 ) sin(tb) (20) 
with the coefficients 7 

=pbi ih xs)(lfl-s)+pal • (\ a x s)(1* • s) 
ci = (paj -ss- pbi)-(lfl x4)-(pal 4- pbi) • (\ b x s)(l a • s) 


and positions in the cube of size 15, centered at the origin, but 
outside of a cube of size 5 centered at the origin, to avoid ini¬ 
tial interference with A. A has 80 vertices, 156 triangles and 234 
edges. B has 85 vertices, 166 triangles and 249 edges. There¬ 
fore, there are 13,280 vertex/triangle, 13,260 triangle/vertex and 
58,266 edge/edge pairs. If we treat the triangle/vertex case as ver¬ 
tex/triangle, there are 26,540 triangle/vertex pairs. 


All types of rejections 



V/T and E/E rejections 



Figure 11: Left: with all types of rejections when objects may 
or mat not collide. Right: with vertex/triangle and edge/edge 
rejections when objects collide. 


c 2 = - pai -(4 x (4 x s))+ pbi -(4 x (4 x s)) 
c 3 =ds-l a *4 

C4 = d (4 x s) • (4 x s) 

( 21 ) 

To find the solution, we perform Newton iterations with initial 
guesses that uniformly sample the search interval at a frequency 
five times larger than b/2n, the period in f(t). Note that since 
be [—71,7t] for pose interpolating screw motion, we may need up 
to three initial guesses. Once a solution is found, a and /3 can be 
computed as 

((bi-q a )x4) • (f(l a ,tb)xljj 

10(4,^) x 41 2 (22) 

q a +a<j>(4,*4)-bi) • 4 

U 2 

We are only interested in solutions where t,a, ft e [0,1]. Note that 
singular situations such as zero length edge or parallel edge may 
be ignored since they are detected by vertex/triangle collision. Fig. 
10 shows an example of a collision between two edges. In the left 
figure, red pluses are the initial guesses and red hexagons are the 
solutions found. Among these solutions, only the ones that satisfy 
0 < < 1 are accepted as valid contacts. Finally, the one with 

smallest t is reported as the collision time. 


The left plot in Fig. 11 shows the timings for all colliding and 
non-colliding cases including the sphere/cylinder rejection. This 
results show that 50% of the cases were rejected by our cylin¬ 
der/sphere object-pair rejection test. The execution time of these 
cylinder/sphere test is negligible. Right plot in Fig. 11 is results 
only for the cases when a collision actually occurred between the 
two objects. In these cases, 50.5% of the vertex/triangle and 66% 
of the edge/edge tests were rejected thanks to the corresponding 
rejection tests discussed earlier. The early rejection results in a 
speedup of 55%. The average execution time for rejection tests 
on 26,540 veltex/triangle and 58,266 edge/edge pairs was 32.8ms, 
which yields about 0.39/is for a single pair. If all the rejection tests 
failed, Newton iteration is performed. The execution time for com¬ 
puting an exact contact points by root finding is 9.5 /is in average. 

Conclusions 

The closed form expression of the trajectory of a point under a 
screw motion is used to predict the collision time and contact point 
between two polyhedra, whose relative motion is approximated 
by one or several screw motion segments. Univariate equations 
for collision conditions are derived for the vertex/triangle and the 
edge/edge collision cases. Several trivial rejection cases are pro¬ 
posed that exploit the properties of screw motions. Our tests show 
that early rejection tests speed up the computation by 55% in cases 
where a collision occurs. 


Results 

To report the execution times and the benefits of our rejection tests, 
50,000 random cases were generated and tested on an Pentium4 
2.5GHz PC. The two polyhedra shown in Fig. 1 were used. The 
object A is placed at the origin at t = 0 and moves along a fixed 
screw trajectory with parameters p = [3,3,3], d = 6, b = 271, and 
s = [1,1, 1 ]/a/ 3. The object B is placed in arbitrary orientations 

6 u • (v x w) = (u x v) • w , u x (v x w) = (u • w)v — (u • v)w , and 
(u x v) • (u x w) = v • w — (u • v) (u • w) 

7 Similarly to (10), these coefficient can be computed after the rel¬ 
evant transformation that aligns s with z axis and ai onto x axis. 
Note that q a and (j)(\ a ,tb) are ai and \ a at time t. 
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